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1.  Introduction 


In  ballistic  performance  testing,  many  statistics  have  been  proposed  as  measures  of  precision. 
Some  of  these  statistics  only  consider  the  spread  of  the  data  in  one  direction,  such  as  the 
horizontal  standard  deviation  and  the  extreme  vertical  dispersion.  Other  measures  attempt  to 
capture  a  sense  of  the  spread  in  both  directions,  e.g.,  extreme  spread  and  circular  probable  error 
(CEP).  Grubbs1  describes  each  of  these  measures,  along  with  several  others  in  his  seminal  self- 
published  reference,  informally  known  as  the  “Red  Book.” 

In  this  report,  we  examine  the  CEP,  defined  here  as  the  radius  of  a  circle  centered  at  the  aimpoint 
which,  in  the  long  run,  is  impacted  by  50%  of  the  projectiles.  We  assume  that  impacts  land  on  a 
two-dimensional  plane  and  that  the  horizontal  (x)  and  vertical  (y)  locations  are  independent.* 
Furthermore,  we  assume  that  the  impacts  follow  a  Gaussian  (normal)  distribution  centered  at  the 
aimpoint.  We  will  first  review  the  cases  in  which  the  only  source  of  variation  in  each  dimension 
is  the  round-to-round  variation,  ox  and  ay.  Then,  we  augment  Grubbs’  approximate  CEP 
formula2  for  those  cases  in  which  additional  error  sources  contribute  to  the  overall  precision  of 
the  weapon  system.  Finally,  we  evaluate  the  accuracy  of  the  approximate  CEP  formula  using 
Monte-Carlo  simulation. 


2.  CEP  for  Bivariate  Independent  Normal  Impact  Data 


2.1  Case  1:  Circular  Data  ( a  =  a x  =  oy ),  Single  Error  Source 

Because  the  impact  locations  are  assumed  to  be  centered  at  the  aimpoint,  we  may  set  the 
aimpoint  to  be  (0,0)  without  loss  of  generality.  Then,  assume  that  X~N( 0,  a2)  and  Y~N( 0,  a2) 

X  Y  X 2  Y ^ 

are  independent  normal  random  variables;  therefore,  -  ~N( 0,1)  and  -  ~  IV  (0,1),  and  —  +  — 

follows  a  Rayleigh  distribution  with  parameter  value  1.  By  definition,  the  CEP  satisfies  the 
following: 

o.5  =  p(V^T^<csp)  =  p(JTE<=)  =  fb(=),  (i: 


'Grubbs,  F.  E.  Statistical  Measures  of  Accuracy  for  Rifleman  and  Missile  Engineers ;  Havre  de  Grace,  MD.  1964. 

In  referring  to  the  two  directions  as  “horizontal”  and  “vertical,”  we  imply  that  the  target  plane  is  upright,  as  is  the  case  for 
many  weapon  systems,  including  calibers  as  large  as  tank  ammunition.  However,  for  artillery  systems,  the  ground  is  the  most 
likely  target  plane,  and  we  instead  refer  to  the  x  and  y  dimensions  as  “deflection”  and  “range,”  respectively. 

2 

"Grubbs,  F.  E.  Approximate  Circular  and  Noncircular  Offset  Probabilities  of  Hitting.  Operations  Research  1964,  12,  51-62. 
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where  FR  is  the  Rayleigh  cumulative  distribution  function.  Since  FR (x)  —  1  —  e  we  have 


0.5  =  1  -  e 


-m/* 


,  which  can  be  algebraically  rearranged  to  obtain  the  following  solution: 


CEP  =  ^J2ln(2)a  *  1.1774er.  (2) 

2.2  Case  2:  Noncircular  Data  ( ax  *  aY ),  Single  Error  Source 

Assume  that  X~N(0,  ox)  and  Y ~/V (0,  a2 )  are  independent  normal  random  variables;  therefore, 
x  Y  x ^ 

- N  (0,1)  and - N  (0,1),  and  Ux  —  —  and  UY  —  —  are  independent  chi-square  random 

ax  &Y  ax  aY 

2 

variables,  each  with  one  degree  of  freedom.  Following  a  similar  argument  given  by  Grubbs,-  we 
start  by  rewriting  the  probability  statement  as 

P(ylx2  +  Y2  <  CEP )  =  P(X2  +  F2  <  CEP2)  =  P(a£Ux  +  a$UY  <  CEP2)  .  (3) 

The  left  side  of  this  last  inequality  is  a  weighted  sum  of  independent  chi-square  random 
variables,  which,  as  Patnaik3  argues,  has  a  distribution  that  can  be  approximated  by  a  chi-squared 
random  variable  when  appropriately  scaled  by  a  positive  constant  a.  This  random  variable  is 
denoted  as  Vn,  where  n  is  the  degrees  of  freedom.  The  approximation  is  found  by  matching  the 
first  two  moments  of  a(axUx  +  <t2Uy)  and  Vn.  That  is,  we  seek  a  and  n  such  that 


E(a(a2Ux  +  a2UY))  =  E(Vn), 

(4) 

and 

Var(a(oxUx  +  aYUY ))  =  Var(Vn). 

(5) 

The  left  side  of  equation  4 

is 

E{a(axUx 

+  oYUY))  =  a(oxE{Ux)  +  aYE(UY))  =  a(a$  +  aY), 

(6) 

while  the  right  side  is 

E(Vn)  =  n. 

(7) 

The  left  side  of  equation  5 

is 

Var(a(<jxUx  +  aYUY ))  =  a2(^axVar(Ux)  +  <jYVar(UY)')  =  2  a2(ax  +  aY), 

(8) 

and  the  right  side  is 

Far(14)  =  2  n. 

(9) 

3 


Patnaik,  P.  B.  The  Use  of  the  Range  as  an  Estimator  of  Variance  in  Statistical  Tests.  Biometrika  1950,  37,  78-87. 
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Since  the  median  of  the  standard  normal  distribution  is  0,  the  numerator  of  the  right  side  of  the 
inequality  needs  to  be  set  to  0  to  find  a  solution  to  the  CEP.  That  is, 


or 


Substituting  for  n,  we  arrive  at  the  final  solution  of 


CEP  = 


(16) 


(17) 


(18) 


2.3  Case  3:  Circular  Data,  Multiple  Error  Sources 

When  multiple  sources  contribute  to  overall  weapon  system  accuracy  (e.g.,  occasion-to-occasion, 
location,  crosswind),  equation  18  still  applies  as  long  as  these  sources  are  uncorrelated  and  the 
vertical  and  horizontal  impacts  are  independent.  However,  in  lieu  of  the  single  error  source,  er, 
we  use  the  total  system  error  calculated  from  the  usual  root-sum-square  formula, 

(JtOTAL  =  V°l2  +02  +  +  °m- 

2.4  Case  4:  Noncircular  Data,  Multiple  Error  Sources 

When  the  weapon  system  accuracy  differs  in  the  vertical  and  horizontal  directions,  the  presence 
of  multiple  error  sources  does  not  change  equation  18  other  than  ox  and  oY  are  replaced  by 
^x— total  =  V^Ji  +(TX2  +  +  and  gy -TOTAL  =  V^ri  +  °Y2  +  "'  +  ffL  respectively. 


3.  Accuracy  of  the  CEP  Approximation  for  Noncircular  Data 


The  accuracy  of  the  approximate  formula  for  CEP  (equation  18)  can  be  determined  using  Monte- 
Carlo  simulation  by  randomly  generating  impact  data  under  various  assumed  values  of  ox  and 
aY*  The  frequency  with  which  the  impacts  fall  within  a  circle  of  radius  CEP  should  be  close  to 
0.5.  The  MATLAET  code  used  to  conduct  this  Monte-Carlo  study  is  given  in  the  appendix. 
Figure  1  is  a  plot  of  values  of  log10(amax/ omin)  vs.  the  frequency  of  hitting  a  circle  of  radius 
CEP,  based  on  100,000,000  simulated  impacts. 


*  . 

If  both  ox  and  ay  are  increased  by  a  common  factor  “k,”  then  it  can  be  shown  that  the  CEP  will  also  increase  by  a  factor  of  k. 
Therefore,  it  is  sufficient  to  just  consider  the  ratio  ax/aY.  In  the  simulation  study,  we  fixed  aY  =  1  and  varied  ax. 

'  MATLAB  is  a  trademark  of  The  MathWorks,  Inc. 
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Since  the  estimated  hit  probabilities  are  all  slightly  higher  than  0.5,  we  see  that  Grubbs’  CEP 
approximation  is  a  slight  overestimate.  The  approximation  has  an  error  between  0.47%  and 
2.17%  when  the  horizontal  and  vertical  errors  are  within  two  orders  of  magnitude  of  each  other. 
The  relationship  appears  to  have  an  asymptotic  value  of  about  0.5073.  A  similar  Monte-Carlo 
study  in  which  the  order  of  magnitude  of  the  error  ratio  was  increased  to  as  large  as  6  upheld  this 
asymptotic  behavior. 


Figure  1.  Relationship  between  relative  size  of  horizontal  and  vertical  errors  and  the  hit  probability  of  a  circle 
centered  at  the  aimpoint  and  having  radius  given  by  Grubbs’  approximate  formula. 
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Intentionally  lelt  blank. 
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Appendix.  MATLAB  Code  for  Evaluating  the  Accuracy  of  Grubbs’ 
Approximate  Circular  Probable  Error  (CEP)  Formula 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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%  OBTAIN  P (HIT)  ESTIMATES  FOR  HORIZONTAL  ERRORS  OF  1  TO  100 


%  KEEP  VERTICAL  ERROR  EQUAL  TO  1 

i= [ 0 : .01:2] ; 
for  j =1 : length ( i ) 
sigx ( j ) =10 Ai  (j) 
sig=sqrt (sigx ( j ) *sigx  ( j ) +1) ; 
v=2* ( ( (sigx ( j ) A4) +1) / (sigA4) ) ; 
CEP=sig* (l-v/9) A1 . 5; 

q, 

o 


%  exponent  values 

%  horizontal  error 
%  root  sum  squares;  NOTE  sigy=l 
%  df  from  chi-square  in  Eq  (6) 


clear  dat 

for  k=l:100  %  generate  100  Ph  estimates 

imp=mvnrnd ( [ 0  0  ]  ,  [ sigx ( 1 ) A2  0;0  1],  1000000);  %  1M  simulated  impacts 

ph ( k) =mean ( imp ( : , 1 ) . A2+imp ( : , 2 ) . A2<=CEPA2 ) ;  %  hit  frequency  on  CEP 

end 

phit ( j ) =mean (ph) ;  %  avg  100  hit  frequencies  to  "incr  n" 

end 


%  GRAPHICS 

plot(i,  phit, ' r- ' ) ;  %  error  exponent  versus  Ph  estimate 

ylabel (' Estimated  Hit  Probability  within  CEP') 
xlabel ( ' log_{ 10 } (\sigma_{max}  /  \sigma_{min} ) ' ) 
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